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We study dynamical response of a nonlinear dimer to a symmetrically injected monochromatic 
wave. We find a domain in the space of frequency and amplitude of the injected wave where all 
stationary solutions are unstable. In this domain scattered waves carry multiple harmonics with 
equidistantly spaced frequencies (frequency comb effect). The instability is related to a symmetry 
protected bound state in the continuum whose response is singular as the amplitude of the injected 
wave tends to zero. 

PACS numbers: 42.65.Sf,42.65.Ky,03.65.Nk,42.65.Hw, 


I. INTRODUCTION 

A nonlinear quantum dimer represents the simplest realization of the discrete nonlinear Schrodinger equation and 
has attracted much interest for a long time [lHlOf. The interest is related to the phenomenon of symmetry breaking 
(self-trapping) ms, mm- On the other hand, non-trivial time-dependent solutions were found in the nonlinear 
dimer [3-0, S]. 

The observation of these remarkable properties of the closed nonlinear dimer implies application of a probing wave 
that opens the dimer. Respectively, temporal equations describing the open nonlinear dimer become non-integrable 
which constitutes the main difference between the closed and open nonlinear dimers. As dependent on the way of 
opening the transmission through nonlinear dimer was studied in Refs. fl2l - l22| where the phenomenon of symmetry 
breaking was reported. What is interesting there is a domain in the space of frequency and amplitude of injected 
wave where stable stationary solutions of the temporal equations do not exist [lj| [20|. Thus one can expect that 
the dynamical response of the nonlinear dimer will display features which can not be described by the stationary 
scattering theory. In particular injection of a monochromatic symmetric wave into the nonlinear plaquette gives rise 
to emission of anti-symmetric satellite waves with frequencies different from the frequency of the incident wave [20| . 
This phenomenon is known as frequency comb (FC) generation and is widely studied in various linear and nonlinear 
systems |23|-l25|. The FC generation can be interpreted as a modulational instability of the continuous-wave pump 
mode (261] revealed even in a single off-cliannel nonlinear cavity [27l |. In this paper we show that illumination of the 
nonlinear dimer with a monochromatic wave in the domain of instability gives rise to FC generation. 


II. COUPLED MODE THEORY EQUATIONS 

The nonlinear dimer can be realized in 2D photonic crystal in the form of coupled microcavities based on two 
identical dielectric defect rods with the Kerr effect [l5| as shown in Fig. |T](a). Taking the radii of the rods small enough 
we can present each rod by a single site variable Aj,j = 1,2 disregarding space inhonrogeneity of electromagnetic 
field in the rods. In terms of the eigenfunctions of the Maxwell equations for each microcavity it means that only the 
monopole mode with the eigenfrequency uiq resides in the photonic crystal waveguide propagation band and thereby 
is relevant in the scattering. For simplicity we disregard the dispersion properties of the waveguide and write for 
the dimer illuminated by light with the amplitude Ei n and frequency to the following temporal coupled mode theory 
(CMT) equations O Him 


—iAi — (ujq + A|Ai|“)Ai -(- 11 A 2 + + A 2 ) — i^jEi n e lut , 

—iA 2 = (wo + A|A 2 1^)A 2 + uA\ -j- * 7 (Ai + A 2 ) — i y /^Ej jn e lult . (1) 

Here the terms X\Aj\ 2 Aj,j = 1,2 account for the Kerr effect of each microcavity, the term is responsible for the 
coupling of the off-channel cavity with the waveguide. The monopole mode of each cavity is localized within a few 
lattice units [30]. If the cavities are positioned far from waveguide we can neglect direct coupling u between them. 
However even for this case the open nonlinear dimer remains cardinally different from the case of the closed dimer 
because of the interaction between the cavities via the continuum. For simplicity we consider the case u = 0. The 
amplitude of the transmitted wave is given by the following equation [lsl 0] 
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FIG. 1: (a) Two identical microcavities made from a Kerr media marked by filled circles are inserted into the square lattice 
photonic crystal of dielectric rods. The ID waveguide is formed by extraction of a linear chain, (b) Two nonlinear sites marked 
by filled circles are positioned symmetrically relative to waveguide and form open nonlinear dimer. 


where Aj(t) = a,j(t)e lult . The open dimer governed by the CMT equations © is shown in Fig. Q] (b). 

Assuming that the solution is stationary a jit) = ajo = const the CMT equations ([1]) are simplified to 

(y + A|ai 0 | 2 )aio + * 7(010 + a 2 o) - i^E in = 0 

{v + A|a 20 | 2 )a 20 + * 7(010 + 020 ) - i^lE in = 0 (3) 

where v = uiq — w. As the input amplitude Ei n increases the system bifurcates from the symmetry preserving solution 
«io = 020 to the symmetry breaking solution I\ ^ / 2 , 6\ — 0 2 = 0,7r where ayo = exp= 1,2 similar to the 
closed nonlinear dimer SSI- However in contrast to the closed nonlinear dimer the open dimer can also transit 
into the phase symmetry breaking solution with I\ = J 2 but 9\ — 0 2 ^ 0, 7r ill 2d]. 


III. INSTABILITY OF STATIONARY SOLUTIONS 


Numerical analysis of stability of the stationary solutions revealed a domain in the space of parameters u> and Ei n 
where all stationary solutions are unstable [15| . Similar result was found in the open plaquette of four nonlinear cites 


10[. In this section we find the domain of instability of stationary solutions of temporal equations (jX I 
To examine the stability of solutions of Eq. © we apply a standard small perturbation technique 

(^) = djo H - {ptj H - Wj)e^ , j = 1>2, 


( 4 ) 


where the second term in Eq. (U) is considered to be small. For the symmetry preserving solutions we have from Eq. 

© 


where according to Eq. © 


aio — <*20 


A q = 


iy/lE in 

v + A/q + 2iy 


I 0 [(u + XI 0 ) 2 + Aj 2 } = jEl, 


and Iq = |A 0 | 2 . Substituting © into Eq. © we obtain the following system of algebraic equations 


-(/* + 2XRe(A 0 )Im(A 0 ))(x 1 - x 2 ) = (v + A I 0 + 2A Im(A 0 ) 2 )(y 1 - y 2 ), 
[n - 2XRe(A 0 )Im(A 0 ))(y 1 - y 2 ) = (v + XI 0 + 2XRe(A 0 ) 2 )(x 1 - x 2 ), 


( 5 ) 
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FIG. 2: (Color online) Domain in space of physical parameters v = ujo — uj and Ei„ where all stationary solutions are unstable 
for 7 = 0.01 (red solid line) and 7 = 0.04 (blue dash line). Other parameters are u>o = 1, A = 0.01. 


and 


-(/x + 2y + 2XRe(A 0 )Im(A 0 ))(xi + x 2 ) = [y + XI 0 + 2XIm(A Q ) 2 {yi + y 2 ) 

(n + 27 - 2XRe(A Q )Im(A 0 ))(yi + y 2 ) = {v + A I 0 + 2XRe(A 0 ) 2 )(xi + x 2 ). (8) 

For Eq. (0 we obtain the eigenvalues 

/i 2 = — {v + XIq){v + 3A/o). (9) 

The equation fi = 0 defines the boundary where the symmetry preserving family of the stationary solutions becomes 
unstable [H]. From Eq. m we obtain equations v + A I c = 0 and v + 3A I c = 0. Substituting these values of I c into 
Eq. ([6]) we obtain for the boundaries of the domain where stable stationary solutions do not exist 

p 2 471/ 

^in \ > 

Efn = ~^y/ 9 + 7 2 ]- ( 10 ) 

The domains of instability of stationary solutions are shown in Fig. [2] Eq. d8]) does not give contribution into 
the domain of instability. Also numerical analysis of the stability has shown that the symmetry breaking and phase 
symmetry breaking solutions fall into the same domain of stability as the symmetry preserving solution. 


IV. NUMERICAL SOLUTIONS OF NONLINEAR TEMPORAL CMT EQUATIONS 

Substituting Aj(t) = aj(t) exp(iwt) into the temporal CMT equations 0 we obtain 

—id\ = (v + A|ai| 2 )ai + *7(01 + a 2 ) - iy/^E in , 

-id 2 = {v + A|a 2 | 2 )a2 + ry(ai + a 2 ) - i^E in . ( 11 ) 

One can see that the solutions possess a symmetry with the half period time shift corresponding to the permutation 
of the sites 


ai(t + T/2) = a 2 (t), a 2 (t + T/2) = ai(t). (12) 

Indeed, after time shift t t + T/ 2 in the first equation in (fill) we obtain the second equation using Eq. m and 
the periodicity of the solutions. Thus, the system of equations HID is reduced to one temporal equation 

- idj = [y + X\aj{t)\ 2 )aj(t) + ij(aj(t) + a,j(t + T/2)) - i^/E in . 


(13) 
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FIG. 3: (Color online) Time evolution of site amplitudes = 1, 2, real and imaginary parts for Ei n = 0.1, v = —0.001, 7 = 

0.04, A = 0.01. 




FIG. 4: (Color online) The same as in Fig. [3] but for the parameters Ei n = 0.4, v = —0.02, 7 = 0.04, A = 0.01. 


Nevertheless the symmetry m does not allow to solve Eq. (fldl) because of unknown period T which strongly depends 
on the intensity of the injected wave. In Fig. [3] and Fig. 0 we present the results of numerical simulations of Eq. 
m in the domain of unstable stationary solutions which demonstrate the symmetry m- Fig. 0also demonstrates 
rachet effect due to absence of the time reversal symmetry in the open dimer. We chose the parameters listed in 
caption of Fig. [3] and Fig. 0 guided by the data on photonic crystal microcavities from Ref. jl5j . 

In order to compare the results with the closed dimer Ref. jlj we present trajectories projected onto the modulus | cij \ 
and phase difference A 9 between cavities in Fig. 0 (a). Although for a small injected amplitude Ei n the trajectories 
look similar to those shown in Ref. |l| however with the growth of the trajectories become asymmetrical relative 
to A 8 —>■ — A 0. The trajectories projected onto the real and imaginary parts of the amplitudes a,j(t) demonstrate the 
most striking difference between the closed and open nonlinear dimer as shown in Fig. 0 (b). While for the closed 
dimer the trajectories form circles centered at the origin of the coordinate system (they are not shown in Fig. 0 (b)) 
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FIG. 5: (Color online) (a) Trajectories projected onto |a, and phase difference Ad = 61—62 and (b) real and imaginary 
parts of amplitudes aj(t),j = 1,2 for different points in the domain of instability: Ei„ = 0.1,^ = —0.001 (blue dash line) and 
Ei n = 0.4, v = —0.02 (red solid line). Other parameters are u>o = 1,7 = 0.04, A = 0.01. 
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FIG. 6: Fourier transform E ou t(f) of the transmitted wave in Log scale for the parameters listed in Fig. [4] 


the trajectories of the open dimer are shifted relative to the coordinate origin. Phase transformation of the injected 
wave Ei n —» E in e la rotates the trajectories in Fig. [5] (b) by the same angle a. 

It is clear that such a complicated time behavior of the amplitudes will reflect at the transmitted wave according 
to Eq. @. Fig. [HI shows the Fourier transformation of the transmitted wave 

E ou t(t) = J dfE out (f) exp(ift) (14) 

that demonstrates sharp peaks spaced equidistantly, FC comb effect. In what follows we define the interval between 
the peaks of the Furrier transform F out (f) as the FC period Hfc- One can see from Fig. [6] \E out (f)\ ^ \E out (—f)\ 
that is a consequence of the rachet effect as seen from Fig. [4j 
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FIG. 7: (Color online) Time evolution of real parts (a) and imaginary parts (b) of symmetric a a = \{a\ + 02 ) (blue lines) and 
antisymmetric mode a a = |(ai — C12) (red lines) for Ei n = 0 . 1 , v = — 0.001 (solid lines) and Ei„ = 0 . 2 , v = — 0.02 (dash lines). 
Other parameters are cjq = 1 ,7 = 0.04, A = 0.01. 


V. THE ASYMPTOTIC EVALUATION OF THE FREQUENCY COMB PERIOD 


The reason for the cardinal difference between the closed and open nonlinear dimers is the symmetry of the system. 
Let us rewrite Eq. CD in terms of the eigenmodes of the closed linear dimer 

-ia s = {v + A[|a s | 2 + 2|a a | 2 ])a s + Aa 2 a* + 2 ija s - i^E in , 

-ia a = ( 1 / + A[|a a | 2 + 2|a s | 2 ])a a + Aa 2 a*. (15) 

where a SiQ = A(ai ± < 22 ) are the symmetric and antisymmetric eigenmodes of the dimer with the eigenfrequencies 

^s,a — ^ 0 * 

Let us take temporarily the dimer linear. The design of the open dimer (Fig. [D implies that the injected wave can 
probe only the symmetric mode with a Breit-Wigner response a s = i^Ei n /(u + 217 ), while the antisymmetric mode 
remains hidden as seen from the CMT equations (1151) . It oscillates with the frequency v however with the uncertain 
amplitude. That defines the antisymmetric mode a a as a symmetry protected bound state in the continuum irt - Col . 
Returning to the site amplitudes we therefore obtain 


a,- = 


f \/7 L'm 
v + 2 fy 


±ae l 


3 = 1,2 


(16) 


making the time behavior of the site amplitudes of the linear dimer non stationary. This equation constitutes the 
time dependent contribution of the bound state in the continuum established for the stationary case in quantum 
mechanical [37] and photonic crystal systems [Hf. 

The nonlinearity results in two effects. The first obvious result is that the resonance eigenfrequency v + A/o of the 
symmetric mode is shifted proportional to Ef n . That agrees with behavior of the instability domain at small Ei n as 
derived in Section III and shown in Fig. [2] The second effect is more sophisticated. For small Ei n the symmetric 
mode a s is almost constant while oscillations of the antisymmetric mode a a are dominant as shown in Fig. [3 As 
seen from the first equation in Eq. m the antisymmetric mode plays the role of a driving force for the mode a s via 
the the nonlinear term A a 2 a*. Then if the frequency of the mode a a is Q then the symmetric mode oscillates with 
double frequency 212 as it is seen from the numerical solution in Fig. [D Respectively, the transmitted wave carries 
the harmonics with the same frequency 212 in accordance to Eq. ©. 

In order to consider these nonlinear effects in the open nonlinear dimer we use the asymptotic methods by Bogoliubov 
and Mitropolsky [38']. Eq. m can be rewritten as follows 


ia s + (v + 2i"f)a s - i^)E in = £F s (a s ,a a ), 
ia a + va a = £F a (a s ,a a ) 


( 17 ) 







7 


where the parameter A is considered as a small parameter e and functions F sa are polynomial functions of a s , a 
determined by Eq. m Then the solution up to the first order in e can be sought in the form 

a s = So(a, (f>) +esi(a,4>) 

a a = a 0 (a, 4>) + eai(a, 4>) (18) 

as functions of the amplitude a and phase </>. They are given by the following equations 


a = eDi(a) 

4> = v + ef2i(a) 

where v is the frequency of oscillations at e = 0. Substitution of Eqs. m, tm> and relation 

. . d(ls,a . 

0’s,a — O —^-h <?—TTJ— ■ 

oa otp 

into Eq. (1151) gives the following equation at the zeroth order in parameter e 


Sq ((X, (/>) 


iy/lEi n 
v + 2*7 ’ 


ao(a, (f>) = aexp(i<f>), 


where the amplitude a is undefined. 

In the first order in e we obtain the following equations 

= ( l> + 2 *7)si + Oo s o + (N 2 + 2a 2 )s 0 
= Wi + SqOq + (a 2 + 2|so| 2 )ao + {iD\ — afli) exp 


One can expand 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Si (a, 0) = J2n F s,n (a) exp (in</>) 

oi(a, <f>) = ^2 n F a , n (o) exp(in4>) . (23) 

According to Ref. j38j there is an uncertainty in choice of functions Si and ai that allows to exclude, for example, 
the first harmonic contributions iq.i, F„ ] that gives the following equations 

Di(a) = 0, ni(o) =a 2 + 2|s 0 | 2 . (24) 

Then, solutions of (1221) are the following 

si (a, <!>) = - ls f + 2 if go + sg exp(2 i<f>) 

ai{a,(j)) = —^jsgexp (~i<j)). (25) 

This equations show that the symmetric solution consists of even terms n = 0, ±2,... in the expansion (P?3l) while the 
antisymmetric solution consists of the odd terms n = ±1,±3,.... The higher orders in the small parameter holds 
these features. From Eq. we have 

(j) = {v + A (a 2 + 2|so| 2 )) t = fl a t (26) 

which yields the FC period S2 _fc = 2f2 a in the first order in A with the amplitude a remaining undefined. This 
amplitude can be determined by the equation a = 0 in Eq. m if the injected amplitude is taken as a small 
parameter e in the perturbation approach. However that approach is successful only in the fourth order in e resulting in 
cumbersome equations. Therefore we estimate the amplitude a averaging the numerical solution over time: a = ( a a (t )). 
The numerical result shown in Fig. [8] (a) is close to analytical result (l26ll in Fig. [8] (b) when the injected amplitude 
is small. Thus, the period of harmonics generated by the open nonlinear dimer can be effectively controlled by the 
injected amplitude. 

Note, the reason for instability related to the bound state in the continuum preserves for u ^ 0 and what is more 
surprising even for different couplings Jj,j = 1, 2 of the sites with the injected wave. For this case the CMT equations 
m will take the following form [29[ 

—id\ = (v + A|ai | 2 )ai + * 7 iai + iyj 717202 - iy/y[Ei n , 

-ia 2 = {v + A|a 2 | 2 )a 2 + iy/j i 7 2 ai + * 7 2 a 2 - iy/y^E in . 


( 27 ) 









8 



FIG. 8: (Color online) (a) Period of harmonics fl vs amplitude Ei n and frequency v of injected monochromatic wave calculated 
numerically, (b) Difference between numerical data and analytical results given by Eq. (1261) . The parameters of the dimer are 
ojo = 1,7 = 0.01, A = 0.01. Below the domain of instability defined by Eq. G0 is shown by red lines. 


By linear transformation 


ai = \/7i(a s + a a ), a 2 = y^n{a s - a a ) 


( 28 ) 


Eqs. (EH) take the following form 

~ia s = va s + 7 + A[(|a s | 2 + 2|a Q | 2 )a s + a 2 a*] + 7 _A[(|a a | 2 + 2|a s | 2 )a Q + a 2 a*] + 2i(^ + a s + 7 _a a ) - iE in , 

-ia a = va a + 7 — A [(| a, s | 2 + 2|a a | 2 )a s + a 2 a*] + 7 +A[(|a Q | 2 + 2|a s | 2 )a Q + a 2 a* + 7 +a 2 a*] (29) 

where = (71 ± 7 2 )/ 2 . One can see that similar to the former symmetric case 71 = 7 2 the mode a a is coupled with 
the injected wave only through the nonlinear terms. 

VI. SUMMARY AND DISCUSSION 

In this paper we considered one of the simplest nonlinear open system, dimer whose closed counterpart is exactly 
integrable system [ij. The term ’’open” means that a linear waveguide is attached to the dimer to allow probing the 
dynamical properties of the dimer. Even in the case of decoupled (u = 0) nonlinear sites they interact with each 
other through the continuum of the waveguide. In the framework of coupled mode theory we examined the stability 
of stationary solutions of Eq. 0 in the parametric space of frequency and amplitude of the probing wave. We found 
a domain where stable stationary solutions do not exist. First such domains were found in open nonlinear plaquette 
[20j together with the associated effect of frequency comb generation. In the present paper we showed a similar effect 
for scattering of a monochromatic wave by a nonlinear dimer. 

The instability of the open nonlinear dimer is related to a symmetry protected bound state in the continuum. When 
the dimer is linear there were two eigenmodes, symmetric and antisymmetric. The symmetrical design of opening of 
the dimer (see Fig. [lj implies that the injected wave couples only with the symmetric mode while there is no direct 
coupling of the injected wave with the antisymmetric mode as seen from Eqs. ED- However owing to nonlinear 
terms in these equations the antisymmetric mode a a is coupled with injected wave through the symmetric mode a s . 
Therefore the antisymmetric mode emerges in the response in the vicinity of the resonance v = 0. 

Numerical solution of the temporal coupled mode theory equations m demonstrates highly nonlinear behavior of 
the site amplitudes cardinally different from the dynamical behavior of the closed dimer in the instability domain. 
Time dependence of these amplitudes holds many harmonics whose frequencies are equidistantly spaced with the 
interval Cl. This interval which defines the FC period Clpc = 2f2 was computed numerically and evaluated by the use 
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of asymptotic methods jUJ to demonstrate an agreement as shown in Fig. [8] Respectively, the injected wave after 
scattering by the nonlinear dimer acquires these harmonics. The value fi goes down with decreasing of the injected 
amplitude that opens a way of all-optical control of the harmonics. 
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